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ABSTRACT 


A model has been developed to estimate the risk 
to people during storage of spent fuel at the reactor site 
and during transportation of the spent fuel "bundles to the 
reprocessing plant. The model begins with the ca-lculation 
of fission product inventojcy in a typical OARDU reactor. 

The fission products decays with time, and the production 
of their daughters has a marked effect on the decaying 
activity of the total inventory. 

Fault trees are developed to compute the accident 
probability in both the Spent Fuel Storage Fool and in the 
transportation of the spent fuel bundles to the reprocessing 
plant. The exposure to people, both due to beta and gamma 
rays following the dispersion of the radioactive plume under 
different weather condition has been estimated. The risk 
to people through external exposure is presented as the 


result. 
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CHAPTER I 


INTRODUCTION 


Nuclear power plants and tHeir fuel-reprocessing 
facilities generate various types of radioactive wastes 
and the safe disposal of these wastes is an integral part 
of the use of nuclear energy .Radioactive waste handling, 
from nuclear power plant to final dumping involves several 
stages and many personnel, as shown in Rigi lil. These 
stages includes cooling of the spent fuel in the Spent 
Fuel Storage Pool (SFSP), transportation to the reprocessing 
facility, transportation to and final dumping at the dumping 
sit e* 

India has a 200 MWe CANDU unit operating and another 
of same capacity awaiting commissioning at Kota, Rajasthan. 
Two units of 200 MWe each are awaiting commissioning at 
Madras, two units of 230 MWe under construction at Narora, 
and another similar pair of units proposed somewhere in 
Gujarat. It becomes obvious that large amount of activity 
will have to be handled in future. Presently the fuel 
reprocessing is done at Tarapur and there are plans to build 
few more of them near different reactor sites. 

Fuel reprocessing for CANDU reactor spent fuel is 
not undertaken in European countries and Canada as it is 
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found uneconomic. These countries have access to uranium 
enrichment facilities for their needs to run either light 
water reactors or breeder programs. However, the case is 
different for India. India requires plutonium as a substi- 
tute to enriched uranium for the subsequent fast breeder 
reactor phase of our nuclear power program.. Plutonium can 
be obtained from natural uranium fuelled thermal reactors. 

As shown later, in Section 2.2, fissile plutonium is pro- 
duced along with the fission products during the reactor 
operation in fuel bundles. This plutonium can be recovered 
by chemically reprocessing the spent fuel. Thereafter, as 
shown in Pig. 1.1, the depleted fuel with fission products 
is stored for an extended period, and finally disposed off 
in geological formations. Any accident or faulty handling 
at any of the stages involved, as shown in Pig. 1.1, poses 
a problem of radioactive release to environment and exposure 
to workers. There is some existing opinion that it is 
necessary to define an acceptable level of safety or in the 
alternative, an acceptable level of risk. This is amply 
brought forward by the IEEE special issue on Nuclear System 
Reliability and Safety [1]. 

1.1 COHOEPT OP RISK ; 

An acceptable risk to the society depends on several 
factors like time, nature and affluence of society, degree 
of technological progress, beliefs and literacy. To assess 
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the societal risk is therefore to confront a very difficult 
problem. According to Otway et al [2], 'The measurement of 
social values and their reconciliation with technical risk 
estimates through the framework of formal decision-making 
methodologies is defined as risk evaluation^ Therefore the 
best method seemingly possible is to compare the risk posed 
by a new technology with the risk posed by an already accepted 
technology. 

To bring clarity to the term 'risk* , Reactor Safety 
Study [3] defines risk to encompass potential fatalities 
and injuries to people and property. Societal risk from 
accidents is defined as 

RISK ^ . 

X Magnitude |; C °SB| . ^g . S |es j 
Another definition given by Kastenberg et al. [4] is 

RISK = [ Probability of occurrence of event] x 

X [Consequence per event ] (1*1) 

This definition in more detailed form is forwarded by 
Smith et al [5] as, 

RISK = AxBxCxDxE (1.2) 
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where, 

A = prohahility of release sequence, 

B = release quantity 

C = relevant measure(s) of the radioactive material 
characteristics (isotopic composition, size 
distribution, solubility, volatility, etc.) 

D =r a measure of environmental transport path 
efficiency 

E = a. measure of population distribution charac- 
teristics and habits 

The essence of the three definitions is however the 
two terms viz. probability and consequence of a particular 
event. However, the consequence assessment is the toughest 
problem since it usually includes injuries (fatal or non- 
fatal) to people and damage to property. The approaches 
prior to Reactor Safety Study handled the problem by convert- 
ing the personnel injuries to lost man hours and property 
damage to monetary penalties. Reactor Safety Study [3] 
however, selected the following four types of consequences, 

a. Early fatalities 

b. Early illness 

c. Late health affects attributable to the accident 

d. Property damage. 
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1.2 RISK ASSESSIIBKT: 

To estimate 'risk' in its totality, we require infor- 
mation of geographical, geological and demographic details, 
accident data at all stages and biological consequences of 
exposure to radiation. To compute risk it is also necessary 
to know the probability of radioactive material leakage 
following an accident. Unfortunately, for the reactors and 
all its associated stages shown in Pig. 1.1, either due to 
complication involved, or due to use of advance engineering 
systems , with inherent high reliability, coupled with precau- 
tions taken to avoid any accident due to public awareness 
and outcry, there is a serious lack of data which makes 
conventional analysis impossible. Therefore the Event Tree 
and Fault Tree methodology are very useful for risk assess- 
ment. Reactor Safety Study [6] uses the above methodology 
to analyze numerous accident chains. For each accident 
chain, a probability of occurrence is computed and the 
resulting consequence is evaluated. The collection of 
probabilities and consequences for the various accident chains 
gave the required point estimates, from which risk calcula- 
tions can be performed. 

1.3 FAULT TREES; 

Fault tree analysis was introduced by Bell Telephone 
Laboratories in 1961 for performing safety evaluations of 
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launch control systems for the Minuteman Program [7], Fault 
tree analysis uses BOOIEAlT AIG-BBEA to depict component 
failures, and through AKI) and OR comhinatorial logic produce 
system failure. Fault tree analysis is a deductive method 
wherehy the analyst starts with the undesired state of the 
system as defined by the TOP event, and then proceeds to 
examine the system, component by component, for these combina 
tions of component fault states that result in the undesired 
system state. 

The relationship of these events to the TOP event are 
described by the fault tree equation by using Boolean operate 
i.e., effectively adding probabilities on 'OR' gates and 
multiplying probabilities on ’AND' gates. These events are 
then individually treated as TOP events and the process con- 
tinued till we reach events for which data is available. Care 
has to be taken to avoid omissions, repetitions and logical 
discrepancies . 

Some of the advantages of Fault tree analysis are [8] s 

(1) Directing the analyst to ferret out failures in a 
deductive way. 

(2) Providing a graphical aid giving system management 
visibility to those removed from the system design. 

(3) Providing options for qualitative or quantitative 
system reliability analysis. 
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(4) Allowing the analyst to concentrate on the particular 
system at a time. 

(5) Providing the analyst with genuine insight into system 
hehaviour . 

The drawbacks listed are 

(1) Its inability to model failures that cannot be decomposed 
through binary logic e.g. system dynamic failure. 

(2) High cost involved and high time consumption, 

(3) Possibility of logical discrepancy in the fault tree, 
and repetition if constructed manually. 

Por a small tree, fault tree construction may look 
attractive. However, for a big multilevel fault tree, the 
process of finding a value or an expression for the resultant 

e 

probability can be cumbersome. Several computer codes are 
available to analyze fault tree quantitatively. In many of 
these codes the numeric probabilities of failure are computed 
by means of Yesely’s Kinetic Tree Theory [9]. One example is 
WAM/BAM [10], which probabilistically evaluates systems 
modeled with Boolean algebra. Another example is the PPEP/ 
KITT [7] code which computes the failure probabilities asso- 
ciated with the fault tree of a system. The minimal cut 
sets of the system (the smallest sets of components that are 
sufficient to cause system failure) which are required by KITT, 
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are provided by the companion PREP. Minimal cut sets are 
also found by a method and algorithm presented by G-arribba 
et al. [11]. 

Improvement on all the above methods has been done by 
Chamow [9]. His method involves the use of directed graph 
and related matrix methods so as to speedily transform the 
graphics into corresponding matrix. The major benefit is 
because the mathematical solutions are readily performed 
by standard matrix techniques, which can be implemented either 
manually or by computer. Other computer aided methods were 
developed by Lapp and Powers [12], Schneewis [13] and 
Semanderes [14], the last one developing the ELRAPT code for 
efficient logic reduction analysis of fault trees. 

Fuel in a nuclear power plant, is located in the reactor 
core, the spent fuel pool, the refuelling process area, and 
the spent fuel shipping cask. By far the largest amount of 
radioactivity is located in the fuel in the reactor core. The 
much larger amount of radioactivity that resides in the core, 
as opposed to other locations is one of the main reason why 
the bulk of attention in safety of nuclear power plants has 
been directed towards potential accidents involving only the 
core. According to Levine and Vesely [15], the refuelling 
operation and the shipping cask are not important contribu- 
tors to overall risk, since it is exceedingly unlikely that 
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the fuel will melt in these situations. They conclude that 
a potentially large release of radioactivity could involve 
the fuel in the reactor core or in spent -fuel pool. 

For any computation of risk at any of the stages mentioned 
above, the prerequisite is calculation of the fission product 
inventories. Clarke [16] in 1972 gave the computer code FISP 
for comprehensive calculation of fission product inventories. 
The FISP program produces number densities and activities of 
fission products for a type of reactor fuel after a given 
irradiation histoiy. The program considers all the radio- 
active and stable nuclides between ‘^Zn and "^Dy to make 
it very comprehensive* The fission product chains upto six 
isotopes long, and half lives above 1 sec. are included. 

The next stepping stone for consequence assessment was 
again provided by Clarke [17]. The WEERIE program is essen- 
tially for assessing the radiological consequences of airborne 
effluents from nuclear installation. The integrated effluent 
is dispersed down wind and PASQUILL weather categories are 
used* The concentration thus obtained are then used to 
obtain inhalation doses. 

More recent work published along this time is the work 
of Fryer and Kaiser [18] in 1978, TIRIOF 2 and TIRIOIT 4 are 
computer programs for use in nuclear safety studies. The 
later one is the improved version. It includes a systematic 



11 


study of plume rise, several refinements of the meteorological 
model employed, a flexible approach to the relationship of 
dose and consequence and an examination of milk ingestion 
path-way. 

Refinement of this work has been carried out by Strachan 
and Godard [19]. Their model calculates the depletion, by 
dry deposition and washout to the ground, from a plume of 
fission products in the atmosphere, and for determining the 
isotopic composition of resulting gamma and beta activities 
in the ground deposit at any selected point downwind. 

The transportation of radioactive materials requires 
special a ttention because of the unique hazard associated with 
such cargoes. The first major woik was attempted in this 
direction by leimkuhler [20].., In the absence of substantial 
experience, specially for shipment of highly radioactive 
materials, the evaluation of risk is quite arbitrary. Leimkuhler 
attempts to solve the problem by making it a safety vs. economy 
problem. He restates the problem as one of cost minimization, 
where total cost includes both the cost of control and the 
liability for accident losses. 

Generally all accident studies consider radioactive 
doses to an essentially stationary population along the 
transportation routes or surrounding the site of an assumed 
transportation accident that has led to dispersal of 
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radioactivity, Yadigaroglu [21] stresses, however, that the 
dose to the traffic passing the fuel-cask-carrying vehicle 
might also be of importance on busy highways. Similarly, 
for highway accidents resulting in some radioactivity 
release, he considers the radioactive dose to the occupants 
of the vehicles that pass in the proximity of point of 
release. 

Smeed and Jacobs [22] have attempted at the analysis of 
fatalities both in the developed and developing countries 
based on the number of road fatalities, vehicles and popula- 
tion. Jacobs was able to highlight the fact that the less 
developed countries (with lower vehicle ownership) will havo 
the higher fatality rates. This fact does havo a bearing at 
calculating the accident probability in roa.d transport. 

More recently, Swindell [23] in 1978 has brought out 
the importance of packaging and shipping cask design in the 
transportation of radioactive material. The interest on 
packaging and transportation was initially evoked in the 
G-eneva IV Conference, by the work of Ferguson and Salmon [24] 
and Richards [25], The special issue of report on packaging 
and transportation highlights the interest created. Devine [26] 
analyses the Department of Transportation [DOT] Regulation 
whose objectives is the protection of the general public, 
transportation workers and property. It emphasises that 
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approaches towards ensuring safety is to establish rigorous 
packaging standards to eliminate the possibility of release of 
contents under adverse conditions or to limit the possible 
release to that quantity that poses no significant hazard. 

1.4 PRESENT WORE; 

The present work is an effort to estimate risk to public 
following an accident involving spent fuel. During the 
operation of the CAEDU reactor, fission products are generated. 
The spent fuel is then left to 'cool’ in the Spent Fuel 
Storage Pool normally for 120 days but in India the cooling 
period is around 400 days. In Chapter II, the methodology 
involved in- calculating the fission products in spent fuel, 
and the decay of activity with time is presented. Chapter III 
contains description of the fault trees used to compute the 
accident probabilities in both the SPSP, as well as in trans- 
portation of the fuel bundles from the reactor .site to the 
reprocessing plant. Monte-Carlo simulation used to quantify 
the fault trees and the confidence interval estimate of the 
TOP event for both spent fuel storage and tia.nsportation are 
given. Chapter 17 contains description on the dispersion of 
the effluent under different weather conditions, and computa- 
tion of dose due to both gamma and beta rays and final risk. 
Chapter 7 contains the results, the conclusions drawn from 
them and an outline for further work. 



CHilPTER II 


FISSION PRODUCT INVENTORY 


The radiation exposure that the public may receive in 
the event of an accident is governed by the amount of fission 
product inventory at that time. The quantity of fission 
products existing can be acciurately estimated from the operat- 
ing history (power level, time, operating schedules). Fission 
products include many atomic species in varying amounts and 
with different physical chemical and radioactive properties. 

Though neutron 'fluence' or 'flux-time' are expressed 
as a measure of exposure the most popular unit of exposure 
of 'burn-up' used is megawatt days per tonne (MI/D/T) of 
uranium. Another unit frequently used is the atom fraction 
bum up. This unit is defined as. 

Number of fission product atoms formed 

p = 

Number of fuel atoms in reactor 

These two units are related by the expression for uranium 235 
fuel, 

= 9.5 X 10^ X p. 


burn up in MWD/T 
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2^1 OPERATING SCHEMES; 

The simplest operating scheme is the batch type in which 
the total fuel loaded in the reactor is discharged at a time. 
The fuel at the centre of the reactor is exposed to a maximum 
neutron flux due to the cosine flux shape. In this scheme, 
the burnup of the fuel is not uniform, and fuel utilisation 
is poor. 

Other schemes to improve the fuel utilisation are 
(l) Graded irradiation, (2) Steady axial movement of fuel, 

(3) Salt pepper method. The steady axial movement method 
is of particular interest since this scheme is utilised in 
the CANDU reactors. Steady axial movement of fuel is 
aimed at achieving uniform burn-ups of the fuel bundles in 
each channel. Fuel bundles in the neighbouring coolant channels 
are moved in the opposite directions. This scheme is, therefore, 
called the bidirectional fuelling. The rate of fuel bundle 
movement is chosen so as to give uniform burn up of all the 
discharged fuel bundles irrespective of their location in 
the reactor. A CANDU reactor discharges 8 fuel bundles per day. 

The burnup achieved in a scheme having bidirectional 
steady axial movement of fuel is far higher than the batch 
type irradiation scheme. Burnup for a typical Rajasthan 
reactor is 7000 MWD/T though theoretically a burnup of 
10,000 MWD/T is possible, due to plutonium 239 build up, 
and its recycle. 
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2.2 FISSION PRODUCT COMPOSITION OF A TYPICAL CANDU REACTOR; 

For the analysis of the spent fuel composition, diffe- 
rential equations governing the various fission product pairs 
in the fuel element are given below. 

d Nj, (25) 

= ^25 ^25 *^25^ (2.1) 

The solution to the above differential gives the rate of 
growth of the fission product pairs from uranium 235, as 
follows ; 

Nj, (25) = N25 (0) (1-exp (-023 e)/(l + cc^^)) (2.2) 

The plutonium 239 fission contribute the following concen- 
tration, 

Wj, (49) = cr^g/(l + a^g) [ (C^ 9 + C2 (l-exp (-<723 

“ (C^ + C2) (1-exp (-0^g Y ©))/o ’49 Y ] 

(2.3) 

Plutonium 240 is not fissionable. Plutonium 241 yields a concen- 
tration of fission products pairs governed by, 

%(41) = [ (1+ a4 i)][C 6© + Cy(l-exp (-0236)7025 

+ (Cg/^49 Y)(l~e2cp (-04gY©)) 

+ ( 09/(740) (1-exp (“<740 6) 

- (Cg + Cr^ + Cg + Cg)(l-exp (-C742 6))/0^^] 

(2.4) 
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Past fission in uranium 238 will produce fission 
products. The concentration of fission products formed 
in fast fission of uranium -238 is governed by. 



( 28 

) = 

(e-l)/(^28-l) [V25 

Np (25) -f 

■ y49 %(49) 




+ V41 Hj, (41) 

n? 

J 


(2.5) 

Burn -up 

of 

fuel is related to fission product pairs 

^7* 


= 

r 

L"! 

, (25) + ^5,(49) + % 

(41) + M 

ij, (28)1 /[K23 (0) 




+ 1^28 (0)1 



(2.6) 

where the 

cons 

tants occurring in 

equation 

(2.3) and 

(2.4) 

are, 







Y 


1 - 

B49 e (1-p) 



(2.7) 



^28 

(0) ^28 '^*^49 ^ 



(2.8) 

02 


^25 

£ (1-p) 1^25(0) 

c^25/('^49Y 

- ^25^ 

(2.9) 


= 

^^28 

(0) 023 “49^ ^^40 ^ 

(1 + «49 

) 

(2.10) 


r= 

C2 

^49 °^49/(^40 “ ^25^ 

“49 

) 

(2.11) 


== 

C3 , 

^4o/(^49^ ” ^^40^ 

^4 ^^40 “ 

^25) ^(“^49^ 








(2.12) 

°6 

=: 

C3 

<=^40/^41 



(2.13) 

O7 

= 

^4 

<^40/ (^41 “ ^25^ 



(2.14) 

=8 

= 

"5 

^(^^41 “ ^49 



(2.15) 

°9 

rr 

( 

+ C4 + C^) 0-43/(043 - 04^) 


(2.16) 
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Tile fission products considered in the study along with 
their half lives is listed in Table 2.2. The yields from 
uranium - 235, Plutonium - 239, and the daughter fission 
products are listed in Appendix A. 

2.3 COMPUTATIOF DETAILS 2 

For the computation of the fission product, the following 


data of a typical CAEDU reactor (RAPP) are 

used . 

1 . 

Thermal utilization factor. 

f = 

0.964 

2 . 

Resonance escape probability 

P = 

0.9309 

3. 

Fast fission factor 

£ = 

1.027 

4. 

Prompt neutrons for every 

•P 25 = 

P 49 = 

2.064 


neutron absorbed 

1.9845 



!1 

H 

2.223 

5. 

Neutron produced per fission 

'-^25 = 

V49 - 

2.426 


reaction 

2.893 



II 

H 

y 

2.317 

6. 

Mean neutron temperature 

T = 

200° C 

7. 

Average burnup 

MWD/T 

= 7000 

8. 

The average cross section, to 

be used, 

are corrected 


Joy ■cis.e-~..a£..¥est.cott • s methods Westcott. substituted.. .the ■■ - 
‘ detailed ■ values of cross section and n.umerically.,.perfoimed - 
the. -integration to obtairL. • the average for a non l/v' “t 3 rpe 
of cross. ..section,, resulting in. 

■ *0 ' ■ = Oq .(g + rs) 


(2.17) 
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where g and s represent the non l/v and spectrum parameters. 
The g and s parameters are tabulated in Table 2,1. The 
parameter r is a measure of the proportion of epithermal 
neutrons in the reactor spectrum. For CAKDUtype reactors 
r = 0.03 in moderator and r = 0.07 in the fuel. The average 
cross-section is obtained by using the following formula^ 

0 = ^ (fn/2) (Tq/T)^ (2.18) 

Using values from Table 2.1, 

g + rs = 0.9482 for all isotopes of uranium 
g + rs = 1.6760 for plutonium 

Equation 2.23 becomes, 

0 = 0.6518795 o' 


The average cross- sections calculated is listed belows 


Isotope 

02000 
( barns ) 

non l/v 
type 

^a, 20000 
(barns ) 

il 

Q 

O 

Q 


o^ = 582 
0^ = 112 

694 

429 

0.1924 

^236 

7 

7 

4.3 


^78 

2.71 

2.71 

1.66 


p239 

0j = 746 
0^ = 280 

1026 

. 1121 

0.3753 

p240 

u 

295 

295 

322 


p241 

u 

0^ = 1025 
Oc = 375 

1400 

1530 

0.3659 
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A total of 79 fission products were considered, as 
listed in Table 2.2, and tlie activities at discharge from 
reactor as well as the decay of the activity at various 
time intervals was computed. The unit of emotivity used^ 
is curie where one curie is defined as 3.7 x 10^^ dis- 
integrations/sec. The activity of a radioactive sample is 
equal to the number of disintegration per second. 

Activity A = N disinfegration/sec . 

_ Q , «69J„ _. ¥ . disintegration/sec. 

j. 

2 

where Ti is the half life of the particular fission product, 

Tj_ for all the 79 fission product is listed in Table 2.2. 

2 

The computer program prepared to compute the activities 
of the’ fission products at discharge, also computes the 
activities of the products as they 'cool down'. The produc- 
tion of daughters of the fission products is also taken 
into consideration. The general solution for the build-up 
of an amount M of the i-th member of the fission product 
chain, 


^i = -^1^2 *•* ^i-1 ^f ^ 




1 - e 


-t 




(2.19) 
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where .|i is the total removal -rate constant and is given by, 

^i = ^i + ^i + c^i 

where , 

= number of atoms of fissionable material 
in reactor 

2 

= fission cross section, cm 

2 -1 

0 = average neutron flux, cm sec 

= yield, atoms of fission product formed directly 
from fission per atom of fuel fissioned. 

^ = radioactive decay rate constant. 

A wide range for cooling time has been taken to analyse 
the activity changes. The time, as also given in Appendix A, 
are 1 hour, 5 hour, 10 hour, 25 hour, 100 hour, 10 days,- 
25 days, 50 days, 80 days, 120 days, 250 days, 1 year and 
2 year. The activities are plotted as a function of time 
in Figure 2.1. The activity decay curve can be devided into 
four ranges. The classification of the fission products is 
primarily done to ease the calculations and later analysis. 

It becomes easy to handle foin* groups of fission products 
than to handle 79 fission products individually. The 
classes can be named as very rapidly saturating (average 
half life = 7.9481 x lO”"^ yrs.), rapidly saturating • 

(average half life = 0.0161 year), slowly saturating 
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(average half life = 0.2402 year) and very slowly saturat- 
ing (average half life = 19^77 year) fission products. The 
first group contains the first 39 fission products, the 
second group contains the next 23 fission products, the 
third group containing the next 9 fission product, and the 
last group having the last 8 fission products. 



Table 2,1; Westcott g and s parameters. 


Element 

T°C 

g . 

Absorption 

s (y=.07) 

'a ^ 

Fission 

s (y=.07)'2 

U 233 


20 

.9983 

1.292 

629.3 

1.0003 

1.222 570.1 

Q 

o 

il 

578 

60 

.9976 

1.380 

632.4 

1.0007 

1.304 571.7 


5 25 

100 

.9972 

1.462 

635.5 

1 . 0011 

1.380 576.3 

lO 


160 

.9971 

1.576 

640.1 

1 . 0019 

1.474 580.7 



220 

.9975 

1.687 

644.6 

1.0029 

1.587 584.8 



280 

.9984 

1.781 

649.1 

1.0040 

1.682 588.9 



^60 

1.0000 

1.906 

655.1 

1.0058 

1.801 594.2 



450 

1.0024 

2.041 

662.0 

1.0081 

1.929 600.2 



540 

1.0052 

2.172 

668.9 

1.0108 

2.053 605.2 



600 

1.0072 

2.259 

673.6 

1.0128 

2.135 610.2 

U 235 


20 

.9780 

.0497 

670.4 

.9759 

-.0502 561.1 

'^ao ~ 

683.04 

60 

.9689 

.0725 

665.3 

.9665 

-.0324 566.4 


577.01 

100 

.9610 

.0923 

660.8 

.9581 

-.0169 552.2 

10 


160 

.9511 

.1159 

655.2 

.9473 

.0017 546.7 



220 

.9433 

.1323 

650.7 

.9384 

.0144 542.0 



280 

.9373 

.1425 

647.1 

.9312 

.0217 538.2 



360 

.9316 

.1492 

643.4 

.9238 

.0251 534.0 



450 

.9273 

.1528 

640.7 

.9177 

.0249 530.5 



540 

.9244 

.1607 

639.0 

.9133 

.0254 528.0 



600 

.9229 

.1639 

638.2 

.9108 

1.0274 526.7 

Pu 239 

20 

1.0723 

2.960 1316.8 

1.0487 

2.329 899.9 


1029.1 

60 

1.1117 

3.085 

1366. 2 

1.0777 

2.432 926.2 

cy 

742.15 

100 

1.1611 

3.130 

1420.4 

1.1150 

2.471 955.9 

^ 10 


160 

1.2582 

3.046 

1514.4 

1.1898 

2.407 1008.1 



220 

1.3836 

2.778 

1624.1 

1.2880 

2.192 1069.8 



280 

1.5345 

2.335 

1747.4 

1.4071 

1.836 1139.6 



360 

1.7658 

1.540 

1928.2 

1.5910,1.195 1242.8 



450 

2.0505 

0.536 

2148.8 

1.8182 

0.384 1369.3 



540 

2.3419 

-.400 

2381.3 

2.0514 

-.376 1503.1 


^r-ntr^ o Orv'T'r? O r\0 T Cm IT 
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Table 2.2; Fission Product Inventory. 


Fission Half 
Product Life 

Name (Yr.) 

Activi ty 
( Curies) 

Fission 

Product 

Name 

Half 

Life 

(Yr.) 

Activity 

(Curie) 

Se®"? 

0.507x10“^ 

0.139x10® 


0.559x10"^ 

0.694 

Se®5 

0.123x10“^ 

0.314x10'^ 

jl35 

0.764x10“^ 

0.510x10^ 

jl36 

0.272x10“^ 

0.637x10^ 


0.105x10“^ 

0.211x10^ 

o 130 
Sn 

0.494x10"^ 

0.142x10'^ 

Sr91 

0.111x10“^ 

0.250x10^ 


0.114x10"^ 

0.587x10^ 

j93 , 

0.117x10“^ 

0.183x10^ 


0.148x10"^ 

0.832x10^ 


0.128x10“^ 

8.48 

As^9 

0.171x10"^ 

0.115x10^ 

1I30 

0.143x10"^ 

1.22 

Se®l“ 

0.35x10“^ 

0.141x10^ 

pa“9 

0.152x10“^ 

0.287x10^ 


0. 418x10”^ 

0.160x10^ 


0.175x10"^ 

15.666 

Sm«5 

0.456x10"^ 

0.255x10'^ 


0.194x10"^ 

0.194x10^ 

Se®3 

0.475x10"'^ 

0.163x10^ 


0.205x10”^ 

1.8389 

-128 

■ 0.475x10“^ 

0.222x10^ 

jl33 

0.237x10”^ 

0.981x10^ 

firs'* 

0.605x10""^ 

0.536x10^ 

pall® 

0.239x10”^ 

0.168x10^ 

hi9s 

0. 989x10”'^ 

0.848x10^ 


0.262x10“^ 

5.0255 

1I34 

0.998x10"^ 

0.275x10^ 

0 121 

Sn 

0. 313x10”^ 

58.881 

3nl28 

0.103xl0‘"^ 

0.120x10^ 

14 3 

Ce ^ 

0.376x10”^ 

0.965x10^ 

Se®^ 

0.108x10“^ 

0.273x10^ 

Br®^ 

0.409x10“^ 

0.0344 

Eu^58 

0.114x10“^ 

0.618x10^ 


0.410x10”^ 

0.291x10^ 

Ba”3 

0.157x10“^ 

0.260x10^ 


0.441x10"^ 

6.6339 

As"^® 

0.173x10“^ 

0.407x10^ 


0.536x10“^ 

0.310x10^ 

CD 

00 

0.239x10“^ 

0.294x10^ 


0.536x10”^ 

0.0791 

BrSJ 

0.273x10“^ 

0.751x10^ 

Cdll5 

0.605x10"^ 

7.587 

Sr9® 

0.308x10“^ 

0.607x10^ 

Mo^^ 

0.759x10"^ 

0.528x10"^ 


0.342x10“^ 

0.113x10^ 

Te^^^ 

0.878x10"^ 

0 

0.366x10^ 

La^^^ 

0.433x10“^ 

0.923x10^ 

rty-l®® 

0.936x10 

0.0222 

Ilu^05 

0.507x10”^ 

0.625x10"^ 


0.103x10"^ 

0.159x10^ 


contd. 



26 


Fission 
Pro duct 
Name 

Half ’ 
Life 
(Yr. ) 

Activity 
( Curie) 

Fission 

Product 

Name 

Half 

Life 

(Tr.) 

Activity 
( Curie) 

Xel53 

0.01445 

0.504x10"^ 

Cd^^^ 

0.1178 

0.1601x10 


0.0189 

0.6469 

Sr®^ 

0.1585 

0.222x10^ 


0.0208 

57.099 

Zr^ 

0.1780 

0.4296 

jl31 

0.0228 

0.102x10^ 

Te'*'^^ 

0 . 2876 

0.0125 

Snl25 

0.0265 

10.024 

Sn^^^ 

0.5726 

42.769 


0.05041 

555.56 


0.7671 

15.205 


0.0550 

O.lllxio"^ 

Eu“® 

1.01 

0.057 


0.0556 

10.146 

Sb^^^ 

2.001 

2.288 

Eu^56 

0.0421 

9.167 


2.600 

0.0291 

Eb®® 

0.0509 

0.559x10“'^ 


4.00 

0.1545 


0.0904 

0,407x10^ 

Kr®5 

10.6 

0.915x10^ 


0.1015 

12.19x10^ 


28.0 

0.974 

Eul°^ 

0.1087 

0.257x10^ 

Sml51 

50.0 

80.0 

1.401 

0.0508 


Time after 
discharge 

Activity 
( Curie) 

Time after 
discharge 

Activity 
( Curie) 

1 hour 

0.559x10^® 

50 days 

0.5152x10° 

£1 

5 hour 

0.5505x10^ 

80 days 

0.2921x10° 
0.145 X 10^ 

10 hour 

0,1999x10^ 

120 days 

25 hoixr 

0.5869x10® 

280 days 

0.1269x10^ 

50 hour 

0.9848xlo'^ 

1 year 

0.4509x10^ 

100 hour 

0.5185xlo'^ 

2 year 

0.6541x10^ 

10 days 

25 days 

0.1425x10^ 

0.8689x10^ 

5 year 

0.3469x10^ 



CHAPTER III 


FAULT TREE CONSTRUCTION AND QUANTIFICATION 


A systematic methodology for the construction of fault 
tree is presented in this chapter. In order to construct 
fault trees- for any system, information is required both to 
describe the system itself, and the operation of the specific 
components -within the system. It is desired to develop 
component models which are independent of any specific system, 
and which depend only upon a set of inputs from other compo- 
nents within the system. With a group of such component 
models available, the description of the system configura- 
tion simply requires a specification of the interconnections 
of inputs and outputs of various components, along with a 
TOP event in order to begin construction of the tree. 

Figure 3*1 gives the logic symbolism and description 
of events used in the construction of fault trees. Boolean 
algebra used for interconnection of the various inputs to 
the top event or describing functional relationship is given 
in Appendix B. 

3.1 FAULT TREE FOR SPEND FUEL STORAGE POOL: 

The fault tree construction process begins with the 
definition of the .undesirable event. Once the undesirable 
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event is selected for the system, it becomes the TOP event 
on the fault tree. For our case the TOP event is the' failure 
of spent fuel storage facility.’ 

Figure 5.2 shows the tree of potential accidents consi- 
dered for the spentful pool. The accident mode is clearly 
divided into two, first the accidents that can be initiated 
by internal plant failures and second, the external forces 
that can potentially cause accidents. 

Release of radioactivity from stored spent-fuel due to 
internal failures can potentially result from heat imbalances 
causing melting of stored fuel or from mechanical damage to 
the fuel assemblies causing release of gap activity. Heat 
imbalance can result from: 

1) loss of coolant water from the spent-fuel storage- 
pool, 

2) loss of the capacity to remove heat from liie pool 
water, which would lead to boiling away of the pool water. 

5) an increase in the heat generation rate in the pool 

because the configuration had been altered into a critical * 
array, again leading to the boil -off of pool water. 

The external causes resulting in failure of spent 
fuel storage facility will be through failure of SFSP con- 
tainment failure. The external causes considered here are 
due to 1) damage due to seismic activity, 2) damage by tornado 
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and 3) damage due to turbine missiles. 

The 12 basic units of the fault tree of SPSP as shown 
in Figure 3.2 are tabulated in Table 3.1 with their failure 
probabilities .per Reactor-Year. The source of data so 
tabulated in Table 3.1 is taken from WASH-1400, Appendix 3t 
4 and 5 [27,28,29]. 

Table 3.1: Fault Tree for Spent Fuel Storage Pool 


Basic unit Probability of 
of fault Accident mode occurrence per 
tree reactor-year 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Dropping of casks due to crane failure 
Human error 

Make-up water unavailable 

Make-up water need not recognised 

Water lost from pool due to inadvertent 
opening of valves 

Dropping of heavy items (failure of 
crane) 


3 X 10 
10 "^ 


-6 


10 

10 


-5 


10 


-3 


3 X 10 


-.6 


Seismic forces resulting in fuel 
arrangement to become critical 

Turbine missiles 

Tornado 

Damage due to seismic forces 
Liquid waste storage tank rupture 
Waste gas storage, tank rupture 


3 X 10 
^-6 


10 

10 


-9 


3 x 10 
-2 


-8 


10 

10 


-2 
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3.2 COMM CAUSE EAIIURE IN SFSP TREE: 

The successful definition of common cause failures 
is necessary to help ensure that all the important contribut- 
ing accident sequences have been defined and that the proba- 
bilities of occurrence of the sequences adequately predicted. 
The identification of common-cause failure therefore is an 
integral part of construction of fault tree. 

Once the TOP event has been defined and the construc- 
tion of fault trees started common cause failures and 
dependencies- are incorporated into the fault trees and their 
quantification in the following five ways, as was also done 
in WASH-1400, Appendix 2 [7]. 

1. The fault trees were constructed to meet the context 
prescribed by TOP event, the fault trees are now conditional. 

2. The fault-trees identify components that are common 
to multiple , systems appearing in an accident sequence, 

3. Bach fault tree is developed to a detailed component 
level in order to locate single component failures and poten- 
tial common-cause failures deep within a system. 

4. Human failures are explicitly included in fault trees, 
and dependencies between human failures were also included 
in the faxilt-tree quantification. 

5. Evaluations, including bounding and sensitivity studies, 
were performed to determine the possible impacts from 


common-cause failures 
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In the SPSP fault tree, it is possible to identify 
common causes and dependencies that are due to not only 
hardware (crane failure leading to dropping of heavy 
items in the pool) but external cause (Damage due to Seismic 
activity) and human failures. Human failure is included in 
the fault tree wherever the operator interfaced with a com- 
ponent or subsystem and could cause failure. Historically, 
human failures were often not included in fault trees j this 
was particularly true when failt tree was constructed at 
the conceptual design stage of the system, where much infor- 
mation was not available. 

In the SPSP fault tree, the common mode failures ax-e; 

(i) Component 1, dropping of casks resulting in damage to 
pool lining and component 6, dropping of heavy objects to 
cause fuel arrangements to become critical. The common 
cause here is failure of the crane. 

(ii) Component 2, dropping of some, heavy items resulting 
in damage to pool lining and component 5» inadvertent 
opening of values in SPSP, resulting in loss of water from 
the pool. The common cause here is human error. 

(iii) Component 7, fuel arrangement becomes critical due to 
seismic forces and component 10, spent fuel pool containment 
failure due to seismic activity. 
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3.3 FAULT TREE OP TTPICAL CANDU SPENT FUEL TRANSPORTATION 
ACCIDENT MODEL: 

A single consignment for shipment in a truck consists 
of 250 fuel bundles out of the 3000 fuel bundles in core. 

So a total of 12 shipment have to be done to transport one 
core load of the fuel bundles. 

For each shipment, the undesirable ’TOP' event is 
'severe/fatal accident of motor vehicle.' Transport Research 
Division, Ministry of Shipping and Transport, Government of 
India in their publication [ 32 ] define severe/fatal accident 
as an accident in which one or more person (occupants of 
motor vehicle) were killed outright or died within 30 days 
as a result of accident. 

Figure 3.3 shows the tree of potential accidents consi- 
dered for the transportation accident model. Like the SFSP 
accident describe in Section 3.1» the tree can be divided 
into two, first the accidents due to inherent failure of some 
critical item in the motor vehicle itself and second, the 
external causes which ultimately result in accident. The 
internal causes considered are defective brakes^ defective 
steering, insufficient or no light, puncture or burst tyre, 
overloading and other serious mechanical defects. The external 
causes resulting in failure considered are defective road 
surface, collision with trains, collision with tree and other 
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stationary objects, and wrong human judgement both on the 
part of driver of motor vehicle and driver of vehicles 
other than the considered motor vehicle. 

The 11 basic units of the fault tree of the transpor- 
tation model as shown in Figure 5*3 are tabifLated in Table 
3.2 with their failure probabilities. The source of data 
tabulated in Table 3.2 are the Ministry of Shipping and 
Transport, Govt, of India publication of Transport Statistics 
[32,33,34]. All the data belong to the basic year of 
1974 - 75 . 


Table 3.2; Fault Tree of Transportation Accident 
Model. 


Basic unit Probability of 

of fault Accident mode occurrence per 

tree kilometer travel 


1 Wrong human judgement, fault of 
driver of motor vehicle 

2 Wrong human judgement, fault of 
driver of vehicles, other than 
the motor vehicle considered 

5 Collision with tree/stationary 

objects 

4 Collision with train 

5 Defective road surface 

6 Insufficient or no light 

7 Puncture or burst tyre 

8 Serious mechanical defects other 

than defective steering or brake 

9 Defective steering 

10 Defective brakes 

11 Overloading 


6.45 X 10 “® 


8.1 

X 

2.33 

X 

2.15 

X 

3.3 

X 

1.2 

X 

2.3 

X 

4,86 

X 

2.5 

X 

5.45 

X 

8.9 

X 


10 ' 


-10 


10 

10 

10 

10 

10 ' 


-8 

-10 

-10 

-10 

-10 


10 

10 

10 

10 


-10 

-10 

-10 

-10 
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3«4 (JEUERAL APPROACH TO PAUIT .TREE QUANTIEICATIOH: 

The fault trees in themselves are qualitative in 
nature, hut they provide a frame work for quantitative proha- 
hilistic analysis. Binary modeling of components in a 
failed state is a characteristic feature of fault tree 
analysis. To facilitate the quantitative analysis it is 
convenient to represent the fault tree in mathematical form, 
and Boolean algebra is an appropriate tool for this purpose. 


Given a fault tree the top event T can he expressed 
aS Boolean functions of primary events E^^ such as 

T = f E 2 , E^, — ) (3.1) 

For fault trees this can he written as, 

T = + M 2 + + (3.2) 

The events are secondary events consisting of intersec- 
tions of primary events, 


m 




where no M^^ is a subset of another V/ith the expression 

in this form, the are termed the critical paths or 
minimal cut sets of the fault tree. 

The probability of T for small probability events 


can be written as, 
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I’(T) = 2^ P (M ) (3.4) 

i=l ^ 

and if primary events are independent then, 

P(M.) = 2. p(0. ) (5.5) 

^ k=l 

In case, common mode failures exist, they are included and 
quantified, as explained in Appendix B. 

3.5 DATA TREATMENT: 

The quantification of fault trees involve two basic 
approaches; point estimation or the random variable evalxia- 
tion. To produce a meaningful point estimate, the input 
data of components should be highly accurate. In actual 
practice, however, extensive failure rate data to back point 
estimation are not usually available. 

The alternative approach is to use the random variable 
technique. Although in use before, this method became very 
popular with the publication of V/ASH-1400 [3]» where it found 
extensive use. 

To express uncertainty in failure rate data a lognormal 
distribution was fitted and the median as well as 5 and 95 
percentile point values are estimated. The lognormal distri~ 
but ion was chosen because of large variability in data, its 
flexibility, its consistency with reliability data and 
because of its property that for a complex system, the 
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pro'ba'bility for the top event is the simple additive of the 
component prohahilities* These justifications are elaborated 
in Appendix B. 

When fault tree is simple, statistical distribution 
algebra can be used to compute the location parameter and 
the variance of TOP event unavailability. When fault tree is 
complicated it is convenient to use 'MOMIE-CARIO' simulation 
of fault tree. 

3.6 POINT ESTIMATION FOB TOP EVENT: 

Using equation 3.5» the top event for the SPSP fault 
tree is denoted by, 

P(T) = X(l) + X(2) + X(3) + X(4) + X(5) + X(6) + X(7) 
+ X(8) + X(9) + X(10) + X(ll) + X(12) + X(13) 

+ X(14) + X(15) (3.6) 

The first 12 probabilities are those listed in Table 3.1, and 
the probabilities X(13), X(14) and X(15) are contribution due 
to commoniE mode failures. G-enerally these are second order 
effects. The computation method for common mode failures, 
which has been dealt in Appsndix B, can be written as. 


X(13) = 

± -.T A 

[X(1).X(6).X(1)]2 = 6 X 10 

(3.7) 

X(14) = 

[X(2).X(5).X(2)]^ = 3 X 10"^ 

(3.8) 

X(15) = 

[X(7).X(10).X(7)]^= 5 X 10"^^ 

(3.9) 
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The point estimation of SPSP comes to be, 

P(T) = 0,02405 (3.10) 

The top event for the transportation model is denoted by the 
Boolean expression, 

P(T) = X(l) + X(2) + X(5) + X(4) + X(5) + X(6) + X(7) 
+ X(8) + X(9) + X(10) + X(ll) (3.12) 

The point estimation of transportation model comes to be, 

P(T) = 0.9249 X 10“^ 

Point estimates for accident probabilities of transportation 
can be obtained by direct evaluation of road fatalities which 
can be established by the use of Smeeds formula [22], 

P/V = 0.0003 (V/P)'"^'^^ 

and Jacobs formula [22] 

P/V = 0.00077 (V/P)”^/^ 

where, 

P = road fatalities, 

V = number of vehicles, 

P = population. 

Jacobs formula is of more interest as it is valid for develop- 
ing countries like India whereas Smeeds formula represents 
developed countries. 
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3.7 MONTE CARLO SIMULATION; 

For each fault appearing in the Boolean expression 
for a system fault tree, a particular failure rate is 
obtained by a random sampling of the appropriate log normal 
distribution. These failure rates are then used to compute 
a value for the top-event failure. This process is repeated 
for a large number of trials to obtain a distribution of top 
event characteristics. Computer code SAMPLE [7] is used for 
Monte-Carlo simulation. 

The SAMPLE program uses Monte-Carlo sampling to 
obtain the mean, standard deviation, probability range and 
distribution for a function, 

Y = f (X^, X 2 , X^) (3.15) 

Given a function Y = f(X^, X 2 , ..., X^), values of the 
location and dispersion parameters of the independent 
variables, and the input distribution, SAMPLE obtains a 
Monte Carlo sampling, X^, X 2 , \ input and 

evaluates the function Y. The sampling is repeated N times 
and the resultant estimates of Y are ordered in ascending 
value y^ y 2 7 ^ "to obtain the limits (percentiles) 

of the distribution of Y. 

The input parameters to the SAICPLE program are N (no. 
of trial runs), median values associated with the 90 percent 
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probability ranges. The upper and lower bounds of the 90 
percent interval are determined by the 95 percent and 5 per- 
cent percentiles of the log normal probabilit 3 ^ distribution. 
The top-event distribution is approximated by a log-normal 
by assigning 90 percent error factor which corresponds to 
the 90 percent interval. The error factor is defined as 
the ratio of the end point to the median value. The program 
prints the result in terms of empirical probability percen- 
tiles. The 50 percent value is the median of the distri- 
bution, and the 5 percent and 95 percent values are bounds 
of 90 percent probability interval. The program also gene- 
rates a frequency function on which the median value and the 
range can be indicated. 



CHA.PTBR IV 


ATMOSPHERIC DISPERSION AND DOSIMETRIC CALCUIATION 

4.1 ATMOSPHERIC DISPERSION: 

As the plume from a site of accident involving 
release of radioactivity is carried away from the site of 
accident by wind, atmospheric diffusion would be continually 
acting to disperse the contaminants at a rate depending on 
the wind speed, thermal stability, and the characteristics 
of the terrain. 

The most important meteorological parameter for this 
analysis is the atmospheric transport vector, or wind. For 
the instantaneous source case, it will determine how quickly 
the effluent puff is carried past*, for the continuous source 
case, the wind will determine the direction of effluent 
travel, and the effluent dilution. Radioactive effluents 
behave essentially the same way as the smoke plume. 

Let *CHI* be the concentration of some effluent as a 
function of space and time. As is proved in Appendix C, the 
effluent concentration 'CHI' in curies/m is Gaussian in 
nature, with standard deviations Oy and Assuming that 

the effluent is emitted at an altitude h from ground level, 
into an infinite atmosphere, the largest value of 'CHI' or 
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the concentration of effluent at ground level ie given by, 

.2 


CHI' = 


Q 


exp ( ) 


TC V Oy 


(4.1) 


2o 


where a , cr^ and h are as described before and v is the 


wind transport velocity. 


Equation (4.1) does not take into account that the 
effluent may be radioactive and may decay while being dispersed. 
This is handled by replacing Q by Q exp (- t) where is 
the decay constant, and t is the time from emission to 
observation. Therefore, equation (4.1) can be written as, 


. -2 

— exp [ - ( t + ^ ) ] 

It V Oy cTg. 2cr^ 

exp [- ( -~ + )] (4.2) 

Tt V Oy Y 2a^ 

The standard deviations Oy and are called the hori- 

zontal and vertical dispersion coefficients and depend on 
atmospheric condition. 

4.2 PASQUILI CONDITIONS; 

To use equations 4.1 and 4.2, Oy and must be 
known as a function of distance from the source. Working 
from experimental data, Pasquill obtained a set of curves 
for o and a for six different atmospheric stability condi- 

y 

tions. They are given in Figures 4.1 and 4.2. A seventh 


'CHI' = 

or 

'CHI’, = 




Distance from source, m 

F ig. 4*1 Horizontal dispersi^ coefficient Oy 
as a functjori of distance from, 

source for the varbus Pasqaitl oandit I 






Hr®E: 

+-U— (44~i-- 


A- Extremely unstable 
B- Moderately unstable 
C- Slightly unstable 
D- Neutral 
E- Slightly stable 
F Moderately stable 


Distance from source , m 


Fig. 4-2 


Vertical dispersion coefficient a 
as a f unc t i on of ‘d istance f rom. ' 
source for the ydrious Pasquill ■ 
ednditiohs. ..'''' ' 
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stability condition, type G, can be approximated by 

V®) = 

The seven stability categories A-G are qiialitatively called 
as the extremely unstable, moderately unstable ^sli^tly 
unstable, neutral, slightly stable, moderately stable, and 
extremely stable categories. The variation of the atmos- 
pheric stability categories during day and night and with 
wind velocity is given below 


Wind speed 
(m/sec) 

Day 

Incoming Solar Radiation 

Right 

Cloud Cover 

Strong 

Moderate 

Slight 

Moderately 

overcast 

Mostly 

clear 

Less than 2 

A 

A-B 

B 

E 

P 

2-3 

A-B 

B 

C 

B 

F 

5-5 

B 

B-C 

C 

D 

E 

5-6 

C 

C-D 

D 

D 

D 

More than 6 

C 

D 

D 

D 

D 


Assuming that all activities involving spent fuel like 
handling and transportation takes place during day-time 
it is clear from above that wind categories to be encoun- 
tered will be A, B, C and D. The effluent concentration has 
been calculated using equation 4«2 for 1, 10, 20 and 50 
kilometers and an effective stack height of 89 metres 
(same as RAPP stack height) for SPSP. Stack height for the 
transportation accident has been taken as 2 metres. 
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4.3 COMPUTATION OP EXPOSURE ANN DOSE: 

Large exposure to radiation for a small time is 
called acute exposure. Continuous exposure to small doses 
is known as chronic exposure. The difference is brought 
out by our model. While the plume is airborne the 'cloud' 
would be dispersed and diluted by turbulent diffusion, I'Haen 
the cloud passed over an area, the people there will be 
exposed to the radiation emanating from the cloud. The 
exposure of people from the passage of the cloud is called 
acute exposure as it would be of relatively short duration 
and the radioactive material would be relatively concentrated. 
Acute exposure is mainly therefore external exposure to 
radioactive plume. 'In- contrast, the exposure of people to 
the radioactive material deposited on the terrain is called 
chronic exposure. The length of exposure here is longer, 
and is mainly through radionuclide ingestion from milk, milk 
products, meat, vegetables and other foodstuff, inhalation 
from resuspended radionuclides may be another path for 
chronic exposure. 

4.3.1 External Exposure to Gama Rays : 

The external exposure due to gamma rays is given 
by following equation [37]: 

1 = 2.62 X 10^ X E X 'CHI' R/Sec. 
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where E is the average energy of all of the gamma -rays per 
disintegration. If ’OHI' is given in Ci/m^ as is for our 
case then, 

i = 0.262 X E z ’CHI' R/Sec. (4.4) 

The derivation of this equation is given in Appendix C. To 
obtain dose equivalent (rems) from exposure (roentgens), X 
must be multiplied by the f-factor and quality factor. The 
f -factor is an energy dependent function which depends on 
the composition of the tissue. For both the factors 

are approximately unity, so, 

H = 0.262 xE X ’CHI' rems/sec. (4.5) 

4.3.2 External Dose From Plume - Beta Rays ; 

As given in Appendix C, the absorbed dose rate to 

air is , 

D . = 4.58 X 10^ X E^ X 'CHI' rad 

air n 

Since the skin is normally exposed to only one-half 
of the cloud, 

■ i) . =r 2,29 X 10^ X Er X 'CHI' rad (4.6) 

air -o 

In equation 4.5, 'CHI’, is in Ci/cm^ , when expressed in 
Ci/m^, i) becomes, 

1) . = 0.229 X Br X 'CHI' 

air -D 


rad 


(4.7) 



CHAPTER Y 


RESUIiTS, COHCIUSIOHS AND SUCGESTIOIS FOR 
FURTHER WORK 

5.1 RESULTS AND DISCUSSION? 

The fission product inventory, and the decay of its 
activity with time is tabulated in Table 2.3. The behaviour 
of the fission products with time is shown by Fig. 2.1, and 
as expected, is found to be in conformity with established 
work [35]. Figure 2.1 shows that the behaviour of the 
fission products, specially due to the production of daughter 
isotopes, changes abruptly at three points, 100 hours, 

120 days and 2 year. This helps to conclude the bouadaries 
for the four groups as described in Section 2.3. The four 
groups mentioned are the very rapidly saturating, lupidly 
saturating, slowly saturating and very slowly saturating groups 

The Monte Carlo simulation was done with the aim 
to arrive at a distribution for the top event. 

5.2 SAMPLE CODE USE ON SFSP TREE; 

The output evaluation for sample size 1200 the dis- 
tribution parameters were found to be 

' Mean = 0.02997 , 

Standard deviation = 0.01306 
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The dose equivaleni; rate in tissue is then given hy^ 

H = 0.229 X 5 X 'CHI’ x f(d, E„ ) rem (4.8) 

B max 

where f is an exjierimentally determined function. The dose 
rate is largest at the surface of the skin, where f = 1> and 
decreases with distance (d) into the tissue. To be on the 
conservative side, we choose f = 1, so, 

6 = 0.229 X Ig X 'CHI’ (4.9) 

The values of ’CHI' calculated as explained in Sec. 4. 2 
are substituted in equations (4.5) and (4.9) together with 
E and Eg and the doses are calculated for the four wind 
condition Jl,B,C, and D at distances of l kilometre, 10 kilometre, 
20 kin.omotre, 50 Kilometer. The total dose due to Gamma 
and Beta rays is tabulated in Appendix D. 

4.4 PINAL RISK ESTIMATION: 

Our objective throughout the study has been to esti- 
mate the risk to population due to accident involving spent 
fuel. The total exposure to population has been estimated 
in rems. The total risk has been estimated by using 
equation (1.1) and (1.2) as described in Chapter I. 

So risk becomes. 

Risk to population = Probability of occurrence of 

accident x release fraction 
X Total dose exposure 


or 


R 


P X f X H 


( 4 . 10 ) 
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H is estimated by using equation (4.5) and (4-9). 

The release fraction f estimation involves the knowledge and 
application of spent fuel can parameters, packing method, 
fraction mechanics, volatility and solubility of fission 
products, and impact analysis. Various experimental studies 
to determine the release fraction has been carried mostly 
in U.S., like dropping of cans from heights of 100 meters 
or more, simulation of impact on computer, actixal accident 
analysis like impact between two trucks, impact against wall. 
But still reliable data about release fraction is yet to te 
established. ¥e solved this problem by calculating the risk 
for 10, 20, 50 and 100 percent release. The risks are 
tabulated in Appendix D. 
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The distribution confidence limits are as follows: 

Confidence 

.lF.ercent) Function 

0.5 

0.01081 

1.0 

0.01197 

2.5 

0.01287 

5.0 

0.014 25 

10.0 

0.01637 

ro 

o 

o 

0.01968 

25.0 

0.02118 

30.0 

0.0221 

40.0 

0.0243 

50. b 

0.027 

60.0 

0.02998 

70.0' 

0.03396 

75.0 

0.03617 

80.0 

0.03861 

90.0 

0.04709 

95.0 

0.05411 

97.0 

0.07195 

99.5 

0.08799 

The frequency 

distribution for the SFSP failure 

plotted in Figure 5. 

1. To analyse the effect of trial : 


on the sample size on the distribution parameters, the sam 
pie size was varied as below; 
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Sample Size 

Mean 

Standard Deviation 

400 

0.030400 

0.01343 

600 

0.030006 

0.01298 

1200 

0.07997 

0.01306 

2200 

0.02993 

0.01286 

3000 

0.02992 

0.01288 


It is concluded from above that the sample size does 
not have much bearing on the distribution parameters. 

5.3 SAMPLE CODE USE ON TYPICAL CAUDU SPENT FUEL 
TRANSPORTATION ACCIDENT TREE; 

The SAMPLE code was run for evaluating the distri- 
bution of the TOP event for this tree. The distribution 
parameters with the sample size is listed below; 

Samnle Size Mean Standard Deviation 


500 

0.2552 

X 

1 

o 

H 

0.4470 X 10"^ 

1200 

0.2374 

X 

10 ^ 

0.3793 X 10“^ 




-6 

—6 

2000 

0.2362 

X 

10 

0.3589 X 10 

3000 

0.2423 

X 

10“^ 

0.36300Z 10"^ 


It is clear that though the best result were obtained 
with sample size 2000, the result was not meaningful. The 
variation of the distribution parameters due to variation 
of the error factors in the input were analysed with the 
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sample size being 2000. As seen from Table 3.2, units 
1 and 3 are the critically important contributors through 
the units 2 and 11 are also important. The effect on the 
distribution parameters with variation of error factors 
is listed in Table 5.1. 


Table 5.1: Effect of error factor on transportation 
tree distribution parameters. 


Component 

1 

Component 

2 

Component 

5 

Compo- Mean 
nent 

4 


standard 

Deviation 

10 

10 

10 

10 

0.2374 

X 

10“^ 

0.3793x10“®' 

20 

. 10 

20 

10 

0.4515 

X 

10~^ 

0.1089x10“^ 

20 

10 

20 

20 

CO 

• 

O 

X 

10"^ 

0.1091x10“^ 

20 

10 

10 

10 

0.5813 

X 

lo-s 

0.9752x10“® 

10 

10 

20 

10 

0 . 2803 

X 

10"^ 

0.6025x10"® 

50 

10 

20 

10 

0.5974 

X 

1 

O 

H 

0.19 34x10“ 5 

50 

10 

50, 

10 

0,6579 

X 

10“^ 

0.2198x10"^ 

5 

10 

5 

10 

0.1517 

X 

10-^ 

0.1312x10“® 

10 

10 

5 

10 

0.2151 

X 

10-® 

0.3545x10“® 

20 

10 

5 

10 

0.3601 

X 

10“^ 

0.966 xl0“® 

5. 

10 

5 

5 

0.1487 

X 

10-s 

0.1304x10“® 


It is clear from Table 5.1» that the best result is 
obtained when the error factor of units 1 and 4 are 5 and 
the rest 9 units are 10. The conclusion drawn here is that 
the probability of components which are critically important 
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The t ransportation accident frequency distribution is 
shown in Pig. 5.2. 

The point estimate of the transportation accident 
probability is 0.9249 x 10 where as using Smeeds formula 
(Eq. 3.14) we arrive at a figure,- 

Probability of accident = 4.68576 x 10“® 
and using Jacobs formula (equation 3.15) we get, 

Probability of accident = 4.3266 x 10~® 

The point estimate is therefore 2.13 times the figure 
given by Jacobs formula which shows that the transporta- 
tion accident fault tree gives a conservative value. 

5.4 EVALUATED RISK RESULTS AND COMCLUSIOUj 

The evaluation of risk to people due to accident 
involving spent fuel, which is the primary objective of the 
study has been evaluated and tabulated in Table 5.2 and 
Appendix L. The risk to people does not reduce exponentially 
with time but the behaviour is similar to the decay of 
fission product inventory as shown in Pig. 2.1. However, 
from Table 5.2, it can be concluded that the risk to people 
decays exponentially with distance. Also it is seen that 
the risk is minimum with wind condition A and progressively 
increases with Pasquill categories B, C and D. We can 
therefore conclude from this, that with unstable wind 
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con'tributors to the TOP event , have to be specified accu- 
rately • With these error factors and sample size 2000 being 
input parameters for the SAJIPIE code, the output distribution 
parameters are: 

Mean = 0,1489 x 10“^ Standard devi. = 0.1304 x 10~^ 

and the distribution confidence limits as given below; 


Confidence (Percent) Function Value 


0.5 

0.2351 

X lO"'^ 

1.0 

0.2612 

X 10“'^ 

2.5 

0.3331 

X 10*’’^ 

5.0 

0.3939 

X 10“^ 

10.0 

0.4905 

X 10"'^ 

20.0 

0.6356 

X lO"'^ 

25.0 

0.6928 

X 10”^ 

30.0 

0.7715 

X 10“^^ 

40.0 

0.9284 

-7 

X 10 ' 

50.0 

0.1109 

X 10“^ 

60.0 

0.1321 

X 10~^ 

70.0 

0.1593 

X 10“^ 

75.0 . 

0.1991 

X lO"^ 



—6 

80.0 

0.2032 

X 10 

90.0 

0.2838 

X 10“^ 

95.0 

0,4015 

X 10“^ 

97.0 

0.4969 

X 10“^ 

99.0 

0.6890 

X 10“^ 

99.5 

0.7939 

X 10“^ 
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Frequency distribution of typical CANDU spent fuel 
transportation accident tree 
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0011(311:1011 the risk to people Is considerably less and 
Increases as the wind condition becomes more stable. The 
risk to population during transportation of the spent fuel 
Is tabulated In Table 5«3. It Is obvious that the magni- 
tude of the risk Is very small and therefore, transportation 
of spent fuel by truck Is quite safe. It can also be 
seen, that the risk to people during transportation of the 
spent fuel by rail (accident probability = 0.14 x 10””^, 
probability of radioactivity release = 0.01 ) Is comparable 
to the risk during transportation by truck. 

5.5 SUGGESTIONS POE PUTURE WORK: 

The present work can be extended to Include: 

1) The list of fission products can be Increased to 
make It more comprehensive. The number of fission products 
considered are 79 but there Is scope for Inclusion of 
other fission products. 

2) The number of decay .chain considered In this work 
was three. But for a better representation of the decay of 
fission product, the number of decay chain should be at 
least 5 or 6. 

5) Sabotage Is an Important component whose Inclusion In 

the fault trees of both SPSP and transportation accident, 
can make the risk analysis more realistic. Usually sabotage 
Is neglected due to lack of adequate data base but the 
Importance Is being recognised now and frequent suggestion 
for Its Inclusion are coming up. 
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4) The effect of terrain on effluent dispersion is 
a major area, where more research effort has to be put. 
This involves more research in the means of geography, 
geology, demography and meteorology. 

5 ) The present work can be extended to include the 
dose due to inhalation of radio nuclides and ingestion 
through foodstuff. 
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Table 5.2 


Risk to people following accident in SPSP. 


Distanc e 

from reactor site = 10 

km. 


Time after 


Weather condition 

discharge 

A 

B 

"C 

D 

1 hour 

0.222x10'“^ 

0.2735 

0.5584x10^ 

0.2168x10^ 

5 hour 

0.3594x10"^ 

0.4429x10“^ 

0.90522 

0.3512x10^ 

10 hour 

0.1307x10“^ 

0.1609x10“^ 

0.32893 

0.1276x10^ 

25 hour 

0.2531xl0~^ 

0.3114x10“^ 

0.6365x10“^ 

0.1234 

50 hour 

0.644191x10“ 

■^0.7927x10“^ 

0.162x10“^ 

0.6286x10“^ 

100 hour 

0.2085x10"'^ 

0.2582x10“^ 

0.524x10“^ 

0.2053x10“^ 

10 days 

0.0350x10“^ 

0.1143x10“^ 

0. 234x10“^ 

0.909x10“^ 

25 days 

0.569x10“^ 

0.700x10""^ 

0.143x10“^ 

0.555x10"^ 

50 days 

0.135x10“^ 

0.3131xl0“‘" 

0.84 r3xl0“^ 

0.927x10“® 

120 days 

0.9488x10“^ 

0.1167x10“'" 

0.2386x10“^ 

0.125x10“® 

280 days 

0.8300x10"'^ 

0,108x10“^ 

0. 2084x10“"^ 

0.809x10“'^ 

400 days 

0.194x10“^ 

0.362x10“ 

0.741x10“® 

O.lSlxlO"''" 

Distance from Reactor Site = 50 km. 


n 

1 hour 

0.9877x10'"^ 

0.395x10"^ 

0. 2236xx 

0.1582x10-^ 

5 hour 

0,1549x10“^ 

0.6399x10“^ 

0.3622x10“-^ 

0.2563 

LO hour 

0.5813x10”^ 

0.232x10“^ 

0.131x10“^ 

o 

0.931x10“^ 

?5 hour 

0.1124x10"'^ 

0.4499x10"'^ 

0. 2547x10“ 

0.1802x10“ 

o 

jO hour 

0.2866x10"^ 

0.114x10“^ 

0.648x10“® 

0.458x10“ 

LOO hour 

0.926x10”^ 

0.370x10“^ 

0 . 209x10"® 

0.148x10“^ 

.0 days 

0.414x10"^ 

0.165x10“^ 

0.938x10“"^ 

0.664x10“® 

13 days 

0. 252x10"^ 

0.93 xl0“^ 

0.572x10“'^’ 

0.405x10”® 

iO days 

0.149x10"^ 

0.596x10“® 

0.337x10“® 

0.239x10“'^ 

.20 days 

0.913x10"’® 

0.113x10“® 

0.913x10"® 

0.873x10“® 

!80 days 

00 days 

0 . 368x10“® 
0.981x10“^ 

0.147x10“^ 

0.524x10“® 

0.335x10“® 

0.996x10“® 

0.391x10“® 

0.812x10“® 
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Table 5.3 

Risk due to transportation of spent fuel by road freight. 


Weather condition 



A 

B 

C C 

D 

Risk at a dis- 
tance from 
accident site 





100 metres 

9.082x10“^° 

1.677x10"^ 

3.082x10”^ 

6.886x10“^ 

200 metres 

1.101x10"^° 

2. 221x10"^° 

3.844x10"^° 

1.006x10"^ 

500 metres 

1.232x10”^^ 

3.539x10"^^ 

7.309x10"^^ 

2.115x10“^° 

1 Kilometre 

7.55 X 10“^^ 

3.906x10“^^ 

1.087x10“^^ 

3.36 xlO“^^ 

2 Kilometre 

3.636x10"^^ 

1,907x10“^'^ 

6.537x10"^ ■ 

1.313x10"^'^ 

Risk due to 

transportation of spent fuel by rail freight. 

100 metre 

1.146x10”^ 

2.117x10"^ 

3.89x10"^ 

8.692x10"’^ 

200 metre 

1.39x10"^° 

2.803xl0~^° 

4.852x10“^^ 

1.269x10“^ 

500 metre 

1.555x10“^^ 

4, 46 7x10”^^ 

9.226x10“^^ 

2.669x10“^^ 

1 Kilometre 

9.5312^0“^^ 

4.93x10"^^ 

1.372x10"^^ 

4,241x10“^^ 

2 Kilometre 

4.59 X 10~^^ 

2.407x10“^^ 

8.252x10"^'^ 

1.657xl0‘'^^ 
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APPENDIX A 


EISSION PRODUCT COMPOSITION 


E . P . Name 

Yield from 
^,235 

(Percent) 

Yield from 

Pu259 

(percent) 

Dau^ter Name 


2 


Br®^ 

Se®5 

1.1 


Br®^ 

ji3e 

5.1 

2.1 

Xe«6 

Sn”0 

2 


Te^^® 

Br®^ 

0.019 


Kr®^ 


0.00035 


ae'’* 


0.056 


33^9 

Sa®l 

0.14 


Br®^ 


0.19 


Ag“^ 

sa^55 

0.033 


Eu155 


0.22 


Br®^ 

jl28 

0.00003 


Xe^^® 


0.92 


Kr®^ 

Mb^® 

0.064 

0.2 

Mo9® 

jl54 

128 

7.8 


Xel54 

c<-u128 

Sn 

0.37 


Sb 

Se®^ 

0.0084 


Br®l 

Bu^58 

0.002 


Ga«® 

Bal39 

6.55 

5.87 

0a^59 

Aa™ 

0.02 


Se^® 

<P 

CD 

<1 

OD 

0.02 


As^® 

Br®^ 

0.51 

0.084 

Kr®5 


5.3 


jS2 


0.011 


O 117 

Sn 


6.4 

5.7 


Eu“5 

0.9 


Ehl°5 

Ga^^ 

0.00011 



ll35 

6.1 

5.7 

10^55 
contd. . . 




I’.P.Nam.e 

Yield from 
jj235 

(Percent) 

Yield from 

■D 239 

Pu 

(Percent) 

Daughter 

Name 


6.3 



SrSl 

5.81 

2.43 

y91 

^93 

6.1 

Zr 

Zr95 


0.0051 



jl50 

0.005 


Xel50 

p^l09 

0.05 

1.4 



0.0078 


Sd^ST 


5.9 

5.5 

Nh^ 

ffai59 ■ 

0.00107 


Eu159 

jl33 

6.6 


0sl55 

PjjllS 

0.01 

0.120 

Cd^-'-^ 

ITb^^ 

0.00061 

0.0036 

Mo®® 

c 121 

Sn 

0.015 

0.045 


r. 145 

Ce 

5.7 

5.3 

Prl43 

Br®^ 

0.00004 


Kr®^ 


0 

3.9 

Ag"°5 


0.0085 


Se^^ 


0.15 

0.37 

Eu153 

Zn^2 

0.00016 

0.00012 

Oa^2 


0.0097 

0.0058 

• lnl^5 


6.06 

6.1 

Tc99 

Pel32 

4.7 

5.1 

Zel52 


0 

0.00006 



0.15 

0.39 

j127 


6.62 

6.69 


P^161 

0.00007 

0.0039 

TV 161 

A 111 

Ag 

0.019 

0.23 

^dUl 

jl31 

3.1 

3.77 

Xe^^^ 


0.0013 

0.071 


m147 

2.7 

2.2 



Contd. . . 



I’.P.Hame 


A-3 


Yield from Yield from Daughter 

^235 Name 

(Percent) (Percent) 



6.35 

5.4 

OeWO 

03^36 

0.0068 

0.11 


Eu156 

0.014 

0.11 

sai57 

Eb®® 

0.00002 

0.00002 

Br®® 


6.0 

5.1 

PI.M1 

Tel®9 

0.35 


jl29 

Eu“5 

5.0 

5.67 



0.0001 

0.003 


Sr®9 

4.79 

1.71 

y89 

Zr^^ 

o 

CM 

• 

5.8 

Nb^ 


0.035 


J.127 

Snl23 

0.0013 


Sbl23 


6.0 

3.79 

EaW4 

Ea^°fi 

0.38 

4.57 


Sbl25 

0.021 


To1®5 


0 

1.94 


Bu«5 

0.033 

, 

Ga^^s 

Kr®5 

0.293 

0.127 

Eb®5 

O 90 

Sr 

5.77 

2.25 

-y90 

137 

Cs 

6.15 

6.63 


o 151 

Sm 

4.4 

0.8 



APPENDIX B 


B.l TOP EVENT AND BOOLEAN ALGEBRA 2 

The TOP event which describes the failure of the 
system, can be expressed as a Boolean algebra function of 
primary events. The following basic properties and laws 
of Boolean algebra are generally used to simplify the 
TOP event function and to eliminate redundant event 

a. Identities ; 

1. A + A = A 

2. A . A = A 

b . Distributive Laws ; 

1 . A(.B + C) = (A.B) + (A.C) 

2. A + (B.C) = (A+B) . (A+C) 

c. Laws of Absorption ; 

1. A + (A . B) = A 

2. A (A . B) = A.B 

Once the TOP event is expressed as a simplified Boolean 
function of primary events, the probability of the former 
can be related to that of later by use of basic operational 
laws to combine probability are given below; 



B-2 


a. UFIOff ; 

Boolean Expression t T = A + 3 

Probability expression 

P(T) = P(A) + P(B) - P(A.B) 

b. INTERSBCTIOE ; 

Boolean expression : T = A.B 

Probability expression, 

P(T) = P(A).P(B/A), if A and B are dependent 

or P(P) = P(A).P(B), if a and B are independent. 

The rare event approximation is applicable when 
the intersection probability P(A.B) is much smaller, than 
the individual probabilities, P(A) and P(B). If A and B 
are independent than P(A.B) = P(A).P(B) and the rare event 
approximation will be valid if cross product from P(a).P(B) 
is much smaller than P(A) and P(B). This will be the case 
when P(A) and P(B) are less than approximately 0.1 [7]. 
Whether the events are independent or dependent and regard- 
less of probabilities, the rare event approximation will 
always give a conservative estimate [5]f whether the 
fault trees are simple or complex, the essential idea is 
to obtain all the minimal cut sets of the TOP event such 
that the TOP event can be written as 



T 


(B.l) 


= M3_ + M 2 + ... 

where the minimal cut-set also termed as critical paths, 
is intersection of primary events 


M. = C-T C. „ C. , . , . C. 

1 il i2 x3 im 


(B.2) 


and no is a subset of another M.. Applying probability 
laws given above to TOP event expression of equation (B.l) 
we get. 


P(T) = 7Z P(Mj_) -II 


hi ^ 

Z P(M. ) P(M.) 
i=l j=l j=l J 


(15^3) 

+ Probability of all possible triple 
combinations 

If rare event approximation is applied then. 


n 


P(T) 


Z. BCM.) 
i=l ^ 


(B.3) 


If primary events are independent. 


m 

P(M.) = Z* B (C.j^) (B.4) 

k=l 

B.2 JUSTIFICATION OF THE CHOICE OF LOG-NORMAL DISTRIBUTION; 

The log-normal is obtained by a transformation of a 
normally distributed variable. It is often used [7] in 
modelling system when factors or percentages characterise 
the input component data variation. If X is a random variable, 
which can vary between X^/f and X^.f, where X^ is some 



mid point reference value and f some factor, then a log- 
normal fits appreciably. 

In this study, the LOG ITORT^AL is suitable because 
the components and the input data in general does vary by 
factors, for example, a potential damage to the SI’S! pool 
due to earth-quake, can have a probability of occurrence 
between 4 x 10 to 2 x 10 per year. Besides giving 
adequate fitting to cases evaluated [7], the log-normal is 
applicable hero from reliability and statistical conside- 
rations. They are; 

1) The positive skewness in the log-normal distribution 
form, can incorix)rate general reliability associated behaviour 
of the data.. The positive skewness accounts for the occurrence 
of less likely but large deviation. 

2) As an a priori distribution the log-normal gives 
coverage to errors which can be skewed towards large values. 

In general, the average value is greater than the most 
probable value, which provides a protective, positive type 
bias, which is retained when the distribution is propagated, 
say as for the calculation of the top event in the fault tree. 

3 ) The log normal distribution being a very flexible one, 
can under the applicable situation, adapt itself as a near 
normal type shape or a near exponential type shape. 



4) The log-normal has an established history of useful 
representation when relative variations (factors) characterise 
the random variable. Its application as a general distribution 
for modelling physical and reliability process is established 
and has been validated. 

5) The assessed data comprises of reliability data in 
the form of probabilities. If the probabilities are decom- 
posed into products of probabilities representing requisites 
for failure, then when the central -limit theorem is applied 
the log normal is the resulting distribution. 

B.3 DEPBNDMCY TREATMENT j 

The handling of dependencies and common mode failures 
while quantifying them is discussed here. As defined earlier, 
a common modp failure is a failure combination (two or more 
failures), which is dependent, or which has more basic 
common causes. In SPSP fault tree the common mode failure 
events are components 1 and 6, 2 and 5 and 7 and 10. The 
following procedure is used to determine the additional 
common-mode contribution of 

Taking the case of component 1 and 6, the contribution 
XQjyj is assessed to be the square root of the product of the 
upper and lower bounds [7, 28] as defined by = X^ 

and Xi^Q^gy ~ ^1 * ^6 X^l values of X(l) 


and X(6). The contribution isfX^^pp^^ %ower 


and hence 
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^GM2 = ^2*^2* V 

^CM3 ~ ^ Zy.Zr^.X^Q 

Since the fault trees include botii-single component and 
common-mode components, it is obvious, that the TOP event 
is dominated by the higher probability, single failure con- 
tributions, and the common mode components are second order 
contributors . 
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C.l EFFLUENT DIFFUSION; 

The diffusion of pollutants is chiefly due to turbu- 
lent eddies in the atmosphere. The individual particles 
become Increasingly separated from one another as the 
result of local atmospheric turbulence and temperature 
variation. The process is called turbulent diffusion. 

Let 'CHI' be the concentration of some effluent as a 
function of space and time. Then 'CHI' is determined by the 
time - dependent diffusion equation, 

[CHI] = (C.l) 

where K is the diffusion coefficient ( cm^/sec. ). Consider 
a pointsource located at the origin of coordinates which 
emits at time t = 0 an instantaneous isotropic puff con- 
taining a total of Q units of effluent into an infinite 
stationary atmosphere. It is shown in pp. 192 [36], that 
the solution to equation (C.l) is 

-r^/4 Kt 

CHI (r, t) = e (C.2) 

(4 n Kt)^/^ 



C-2 


where r is the distance from the origin C.2 is based on the 
assumption that the atmosphere is isotropic. >/hen the atmos- 
phere is non-isotropic , the solution is 

'CHI' (x,y,z,t) = ^ 

(4 (K^Kp^ 

X exp [“ i (f“ + |- + |- )] (C.3) 

X y z 

The movement of effluent in the direction of the wind is 
due to the wind itself and not due to diffusion. So the diffu- 
sion in x-direction is ignored. The physical picture is that 
of a thin disk which is moving along the x-axis with wind 
speed V, which is continuously spreading out in y and z 
direction as it moves. So equation (C.;^) is reduced to 

'CHI' (x,y,z,t) = ::: x“ 

rtvt 

exp [- 4^ ( )] (C.4) 

It can be seen that the effluent is Gaussian ■ in the y- and 
z-direction, with the standard deviations and gi'^en by, 

2a^ = 4tEy , 2^2 = 4tK^ 

So, 



'CHI' (x,y,z) 


0.3 


Q 

2 AV 



)] (C.5) 


It was assumed till here that the effl^e^i'fcs are omitted 
at the origin of coordinates. To consider emission at an 
altitude h, the method of images was applied. So for the 
concentration at ground level, 2=0 


' CHI' 


Q 


exp 



(C.6) 


The value of 'CHI' is largest along the centre line of the 
plume where y = 0. The concentration there is, 


' CHI' 


Q 

XV Oy cr^ 



(C.7) 


C.2 EXTERNAL DOSE FROM PlUI^E; 

C.2.1 Gama -Rays; 

To compute the ffama-ray exposure from a specified 
radiation field, it is necessary to determine the energy 
ahsorhed from the Sama-rays in the air, and then convert it 
to roentgens, where an exposure of corresponds to an 

energy deposition of 

IR = 5.47 X lo'^ MeV/g 



C.4 


Since 1 MeV = 1,60 x lO”"^ ergs, it follows that 

1 R = 5.47 X 10*^ X 1.60 x 10“^ 

= 87,5 ergs/g. 

The total collision density at a point where G-aimna-ray 
intensity is I is given by, 

R = I p 

where p is the total attenuation coefficient. If the gamma- 
ray was absorbed at each collision, then the rate at which 
energy is deposited per unit volume in the medium will be 
BF = El ji , where E is the energy of gamma-rays. 

With both the photo electric effect and pair produc- 
tion, the incident photon is absorbed and most secondary radiu' 
tion - the X-rays, electrons and annihilation radiation is 
also absorbed in the medixm. In Compton scattering, the 
only energy deposited is the kinetic energy of the recoiling 
electron. The energy deposition rate per unit volume by 
Compton scattering is E I where p is the Compton 
absorption coefficient. So the total energy eposition rate 
W per unit volume can be written as 

® I <l"pe + + I'ca* 

= 3 I Pa 


(0.9) 



0.5 


where jj.^ = M-ca called the energy absorp- 

tion coefficient. So it follows from equation (0.8) and 
(0.9) that the exposure rate X is, (where P denotes density of 

X = I B / 5.47 X lo"^ 

0 , 

= 1.85 X 10"® I E R/Seo. 

3 . 

If the gama-ray buildup flux is denoted by then the 

above equation can be rewritten as 

X= 1.85x10“® B(-p^/P)^^^ R/Sec. (C.IO) 

It has been shown in [56] that the build up flux at the 
centre of the cloud is 

s 

^yb 

where S is equal to 5.7 x 10^® 'CHI’ where the radionuclide 
concentration in the cloud is 'CHI ' Ci/cm^ . 

Therefore equation (C.IO) becomes, 

.V . n^-3 5.7 X 10^® X 'CHI' „/a 

X = 1.85 X 10 2^ ^ 

= 2.62 X 10® X 'CHI' X E R/Sec 

where P = 1.295 x 10"® g/cm^ is the density of air. If 

more than one gamma-ray is emitted by the nuclido 

X = 2,62 X 10® X 'CHI' X I R/Soc. 


C.6 


If *CHI’ is given in Ci/m^ then, 

X = 0.262 X 'CHI’ X I (C.ll) 

The dose equivalent can be found by xHtiltiplying 
above equation by f-factor and quality factor. However, 
since both the factors are approximately unity, 

H = 0.262 X 'CHI' X I 

C.2.2 Beta-Rays j 

At any point in an infinite cloud, as many beta-rays 
are absorbed per cm^/sec. as are emitted. If the cloud con- 
tains a radionuclide of 'CHI' Ci/cm^ emitting an average of 
Eg MeV per disintegration, the total amount of energy 
absorbed is 3.7 x 10^^ x'CHI' x Eg x 1.6 x 10~^ ergs/cm^— sec. 
This energy is absorbed in air having a mass of 1.293 x lO”^ 
g/cm^, so that the energy absorption rate is 

3.7 X 10^° X 'CHI' X Eg X 1.6 X 10“^ 

1.293 X 10 ^ 

= 4.58 X lo”^ X 'CHI' X Eg erg/g-sec. 

The absorption of 100 erg/g gives an absorbed dose of 1 rad. 
Thus the absorbed dose rate to air, 

^air ~ x 10^ x 'CHI' x Eg rad/sec. 

Since the skin is exposed to only one-half of a cloud, 



0,7 


i) . = 2.29 X 10^ X ’CHI' X rad/sec. 

air n 

where 'CHI' is expressed in Ci/m^, becomes, 

i) . = 0.229 X 'CHI' X Et, 

air -D 

The dose equivalent rate in tissue is then given by, 

H = 0.229 X 'CHI' X Eg X f 

where f is an experimentally determined function. The dose 
rate is largest at the surface of the skin where f = 1, and 
decreases rapidly inside the tissue. To bo conservative, 
therefore, skin dose is computed with f = 1. 
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DOSE AND RISK TO PEOPLE AT A DISTANCE OP 1 ZILOMETRE 

PROM SPSP 


Time after 
discharge 


Weather 

Condition 


A 

B 

C 

D 

1 hour 





Dose 

O.l^Txlo"^ 

0.625x10^ 

0.865x10'^ 

0.141x10^ 

50 ^ risk 

0.185x10^ 

0.844x10^ 

0.116x10^ 

0.190x10^ 

205 ^ risk 

0.743x10^ 

0.337x10^ 

0.467x10^ 

0.763x10^ 

5 hour 





Dose 

0.222x10^ 

0.101x10^ 

0.140x10^ 

0.229x10"^ 

50 ^ risk 

0.301x10^ 

0.137x10^ 

0.189x10^ 

0,309x10^ 

20^ risk 

0.120x10^ 

0.546x10^ 

0.757x10^ 

0.124x10^ 

10 hour 





Dose 

0.810x10^ 

0.367x10^ 

0.509x10^ 

0.832x10^ 

50 ^ risk 

0.109x10^ 

0.496x10^ 

0.687x10^ 

0.112x10^ 

20^ risk • 

0.4375 

0.198x10^ 

, 0.275x10^ 

0.449x10^ 

25 hour 





Dose 

0.156x10^ 

0.712x10^ 

0.985x10^ 

0.161x10^ 

50?b risk 

0.2116 

0.9612 

0.133x10^ 

0.217x10^ 

20^ risk 

0.847x10“^ 

0.3844 

0 . 5323 , 

0.8498 

50 hour 





Dose 

0.399x10^ 

0,181x10^ 

0.250x10^ 

0.409x10^ 

505'!^ risk 

0.538x10"^ 

0.2446 

0.3387 

0.5534x10“ 

20^ risk 

0.2154x10“^ 

0.978x10"^ 

0.1354 

0. 221x10“^ 

100 hour 





Dose 

0.129x10^ 

0.586x10^ 

0.8116x10^ 

0.1326x10^ 

50^risk 

0.174x10'’^ 

0.791x10"^ 

0.1095 

0.1791 

205 ^ risk 

0.697xl0~^ 

0.316x10“^ 

0.0438 

0.716 

10 days 





Dose 

0.5776 

0.262x10^ 

0.363x10^ 

0.5934x10^ 

50 ^ risk 

0.779x10"^ 

0.354x10'"^ 

0.490x10“^ 

0.801x10"^ 

20^ risk 

0,311x10“^ 

0.141x10^^ 

0.190x10”^ 

0.320x10“^ 




contd. . . 


i 


Time after 
discharge 


Weather Condition 


25 days 

Dose 

50^ risk 
2 % risk 
50 days 

Dose 

50^ risk 
20^ risk 
120 days 
Dose 

5C^ risk 
20^ risk 
280 days 
Dose 

50^ risk 
20^ risk 

400 days 

Dose 

505^ risk 
20;^ risk 

T?T.«=ik to -peoi 

1 hour 

Dose 

5C^ risk 
2C^ risk 
5 hour 

Dose 

5C^ risk 
20^ risk 
10 hour 

Dose 

50^ risk 
20^ risk 


0.3524 0.160x10 

0.475x10“^ 0.216x10”^ 
0.190xl0~^ 0.864xl0~^ 


0.160x10^ 0.221x10 0.36 2x10- 

0.216x10”^ 0.299x10"^ 0.488x10' 


0.907x10“-^ 0.9443 
0.280x10“^ 0.127x10"^ 
0.112x10"^ 0.510x10"^ 


0.587x10"-^ 0.2669 
0.293x10”^ 0.95 X 10“^ 
0.149x10”^ 0.42 X 10”^ 


0.514x10“^ 0,233x10“^ 
0.694x10"'^ 0.315x10“^ 


0.119x10"'^ 0.195x10 ^ 


0.9307 ^ 0.9136 

0.676x10“^ 0.9884x10^ 
0.346x10”^ 0-. 496x10” 


0.3696 0.6039 ^ 

0.498x10"^ 0.8152x10" 
0.199x10"^ 0.293x10” 


0.323x10" 0.528x10 

0.436x10”^ 0.713x10”^ 


0.277x10""^ 0.126x10”^ 0.174x10"^ 0.285x10” 


0.182x10"^ 0.8296x10”^ 
0.246x10”"^ 0.112x10”^ 
0,986x10”^ 0.448x10""^ 


0.114x10”^ 0.1877x10” 
0.155x10”^ 0.253x10”^ 
0.820x10”^ 0,101x10 


lie at a distfl.nce of 20 Km front SFSP 


0 . 2926 


0.395x10' 

0.158xl0‘ 


0.474x10 

0.639x10" 

0.255x10' 


0.172x10 

0.232x10"^ 

0. 930x1 0“4 


0.395x10^ 

0.7465x10^ 

0.357x10 

0.533x10”^ 

0.1007x10^ 

0.483x10 

0.2133x10”^ 

0.4031 

0.193x10 

0.6398 

0.120x10^ 

0.579x10 

0.863x10”^ 

O.I632 

0.7824 

0.345x10"^ 

0.652x10“^ 

0.3129 

0. 2325 

0.439x10^ 

”1 

0. 290x10 

0.313x10”^ 

0.593x10" 

0.284 3 

0.1255x10"^ 

0.237x10" 

0.1137 


contd. . . 



D-3 


Time 


Weather Condition 

discharge 

A 

B 

0 D 

25 hour 

Dose 

0.333x10“^ 

0.449x10“^ 

0.8501 0.4075x10^ 

50 ^ risk 

0.449x10 ^ 

0.607x10“^ 

0.1147x10"^ 0.5502x10"^ 

20^ risk 

0.179x10“^ 

0.242x10“^ 

0.459x10~^ 0.220x10“^ 

50 hour 

Dose 

0.848x10“^ 

O.mxlO"-*- 

0.2163 0.103x10^ 

50 % risk 

0.114x10“'^ 

0.154x10“^ 

0.292x10“^ 0.140x10“^ 

20 fo risk 

0.458x10“^ 

0.618x10“^- 

0.116x10“^ 0.560x10“^ 

100 hour 

Dose 

0.274x10“" 

0.37 xlO“^ 

0.699x10“^ 0.3355 

50 % risk 

0.370x10“^ 

0.50 xlO“^ 

0.944x10“^ 0.457x10“^ 

20 ^ risk 

O.USxlO"*^ 

0, 2 xlO“^ 

0.377x10“^ 0.181x10“^ 

10 days 

Dose 

0 . 122 x 10 "^ 

0 . 165 x 10 “ 

0 . 313 x 10 “-^ 0.1501 

50 fo risk 

0 . 165 x 10 ”^ 

0 . 223 x 10 “^ 

0 . 422 x 10 ”^ 0 . 202 x 10 "'^ 

20 ^ risk 

0 . 663 x 10 “^ 

0 . 895 x 10 “^ 

0 . 169 x 10 “^ 0.810x10"^ 

25 days 

Dose 

0.749x10“'^ 

0 . 101 x 10 “^ 

0 . 191 x 10 “-^ 0.916x10“-^ 

505 ^ risk 

0 . 101 x 10 “^ 

0.136x10““^ 

0 . 258 x 10 “^ 0 . 123 x 10 "^ 

20 fo risk 

0 . 404 x 10 “^ 

0. 546x10“^ 

0 . 103 x 10 “^ 0 , 494 x 10 “^ 

50 days 

Dose 

0.442x10”^ 

0 . 596 x 10 “^ 

0.11 2x10“ 0.540x10"'^ 

50 ^.:risk 

0 . 596 x 10 “^ 

0 . 805 x 10 “^ 

0 . 152 x 10 "'^ 0 . 729 x 10 “^^ 

20 % risk 

0 . 238 x 10 ”^ 

0 . 322 x 10 "^ 

0.608x10“^ 0 . 291 x 10 ”''^ 

120 days 

Dose 

0 . 109 x 10 "^ 

0 . 147 x 10 “'^ 

0.773x10“^ 0.299x10“^ 

50 % risk 

0 . 147 x 10 “'^ 

0.199x10“^ 

0.104x10“^" 0.404x10“^ 

20 ^ risk 

0.59 xlO”® 

0.796x10“'^ 

0.417x10"^ o.ieixio""^ 

400 days 

Dose 

0 . 388 x 10 “° 

0 . 524 x 10 “^ 

0.99 xlO""^ 0 . 474 x 10 ’^ 

50 % risk 

0 . 524 x 10 “® 

0 . 707 x 10 “^ 

0.133x10"^ 0.648x10“^ 

20‘/o risk 

0 . 209 x 10 “® 

0 . 283 x 10 “*^ 

0 . 534 x 10 “^ 0 . 256 x 10 “^ 
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